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ABSTRACT 

Radio-wave scattering is enhanced dramatically for galactic center 
sources in a region with radius ^ 15 arc min. Using scattering from Sgr A* 
and other sources, we show that pulse broadening for pulsars in the Galac- 
tic center is at least 6.3 z/ -4 seconds [y = radio frequency in GHz) and is 
G\ ■ most likely 50-200 times larger because the relevant scattering screen ap- 

pears to be within the Galactic center region itself. Pulsars beyond — but 
viewed through — the Galactic center suffer even greater pulse broadening 
and are angularly broadened by up to ~ 2 arc min. Periodicity searches 
at radio frequencies are likely to find only long period pulsars and, then, 
only if optimized by using frequencies ^ 7 GHz and by testing for small 
numbers of harmonics in the power spectrum. The optimal frequency is 
v ~ 7.3 GHz(Ao.iPv^)~ 1//4 where Aq.i is the distance of the scattering 
region from Sgr A* in units of 0.1 kpc, P is the period (seconds), and a 
is the spectral index. A search for compact sources using aperture syn- 
thesis may be far more successful than searches for periodicities because 
the angular broadening is not so large as to desensitize the survey. We 
estimate that the number of detectable pulsars in the Galactic center may 
range from < 1 to 100, with the larger values resulting from recent, vigor- 
ous starbursts. Such pulsars provide unique opportunities for probing the 
ionized gas, gravitational potential, and stellar population near Sgr A*. 

Subject headings: scattering — pulsars: general — Galaxy: center 
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1. INTRODUCTION 

Only 15 radio pulsars of the current catalog of ap- 
proximately 700 are within 5° of Sgr A 1 and none are 
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likely to find many pulsars. However, some pulsars 
have been detected at frequencies as high as 30 GHz 
(Xilouris et al. 1995; Beiradakis et al. 199!:; [zvekova| 



within 1° (Taylor, Manchester, & Lyne 1993). DifT 
covering pulsars in or near the Galactic center (GC) 
would provide exciting opportunities for probing mag- 
netoionic material, the gravitational potential, and 
star form ation in the GC (e.g., Hartmann 1995] ; So- 



et al. 1994; Kramer et al. 1994; Malofeev et al. 1994) 



|fuc 1994 ). Counterparts to X- and 7-ray GC sources 



may also involve radio pulsars or their progenitors 
Moreover, dynamical constraints on the mass in the 
GC combined with infrared observations suggest that 
neutron stars may abound there. 

It may be impossible to find radio pulsars in the 
GC as sources of periodic radio emission, however. 
The difficulty arises because radio-wave scattering, 
which is particularly severe from sources in the GC 
' e.g., Sgr A *, |Backcr et al. 1993j ; |Lo et al. 1995|: 



Backer 198S ; and several OH/IR stars, Frail et al. 1994; 
van Langevelde et al. 1992), broadens pulsar pulses by 
extraordinary amounts. The actual level of angular 
broadening toward the GC at 1 GHz is ~ 10 times 
greater than that predicted by a recent model for the 
distribution of free electrons in the Galaxy ( Taylor & 
Cordes 1993j ; hereafter TC), even though this model 
includes a general enhancement of scattering toward 
the inner Galaxy. For pulsars in the GC, the resul- 
tant pulse broadening is at least 100 times larger than 
predicted by the TC model. The TC model purposely 
excluded an explicit treatment of the electron density 
in the GC, primarily because it would have been, at 
best, poorly quantified. More recent data on angu- 
lar broadening of GC sources, however, now allows 
us to estimate the amount by which pulse broaden- 
ing should b e augmented. Our results complement 
recent work (Johnston 1994) that suggests a paucity 
of pulsars inside of the molecular ring at Galactocen- 
tric radius ~ 4 kpc. Johnston (1994) considered a 
large-scale region toward the inner Galaxy and did 
not discuss the immediate vicinity of the GC. 

Pulse broadening is so large in the GC that, at 
conventional frequencies used for pulsar searches, 0.4- 
1.4 GHz, the pulsed flux is orders of magnitude too 
small to be detected, as we show below. The strong 
frequency dependence of pulse broadening (oc v~ A ) 
suggests that frequencies at much higher frequencies, 
5-20 GHz, will mitigate the effects of scattering. But 
pulsar spectra generally decline at higher frequen- 
cies, sometimes precipitously, making this option un- 



Though these are relatively nearby objects, their de- 
tections suggest that some pulsars may show sufficient 
spectral flattening that they would be detectable if as 
far as the GC. 

In this paper we discuss pulse broadening from 
sources in and beyond the GC. We show how to 
optimize a periodicity search in terms of radio fre- 
quency and number of harmonics. We also discuss 
searches for pulsars as steady sources rather than pe- 
riodic ones. Namely, we show that aperture synthesis 
surveys using the Very Large Array (VLA), for ex- 
ample, should be successful in finding objects with 
pulsar-like spectra and polarization. 

In we detail the assumed geometry and angular 
broadening. In §^ we derive the resulting pulse broad- 
ening timescales. §^ discusses the pessimistic implica- 
tions of pulse broadening for pulsar surveys involving 
periodicity searches. In ^ we argue that aperture 
synthesis surveys may be more successful than peri- 
odicity searches in identifying radio pulsars in the GC. 
In we estimate the likely number of pulsars in the 
GC. In §[7] we present our conclusions and recommen- 
dations. 

2. ANGULAR BROADENING AND THE 
SCATTERING GEOMETRY 

The angular diameter of Sgr A* is 1.3" at 1 GHz 



and scales as 



in accord with interstellar scatter- 



ing ( Lo et al. 1993 ). OH/IR stars in the vicinity of 
the GC also show very heavy scatt ering, at least for 
those stars within ~ 15' of Sgr A* ([Frail et al. 1994 ; 
van Langevelde et al. 1992 ). We use the measured 



angular broadening to predict the pulse broadening 
expected from pulsars in the GC, but with the follow- 
ing caveat: we must assume the location along the line 
of sight of the scattering "screen" that enhances the 
scattering. As van Langevelde et al. (1992) discuss, 
the location of the screen is by no means determined 
and has implications for the depth of free- free absorp- 
tion to the GC and the outer scale of the electron- 
density variations. For now, we treat the location of 
the screen as a parameter of the overall geometry. We 
constrain its location using recent scattering measure- 
ments for the GC region (Lazio, Cordes & Frail 1996). 
We model the scattering as arising from a discrete 
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screen, as in the top of Fig. [|. Let 9 S be the scat- 
tering angle produced by the screen and 9 Q the ob- 
served scattering angle. In the following we assume 
the small-angle approximation: 8 s ,0 o <C 1. For a 
source at infinite distance, 9 a = 9 S because the waves 
incident on the screen are plane waves. Sources at fi- 
nite distance are scattered less due to wave sphericity. 
If the source is at distance D and the screen is in front 
of the source by distance A, the observed scattering 
angle is 9 D — (A/D)9 S . If the scattering screen is 
very close to the GC, then A -C D and the observed 
scattering diameter of an extragalactic source is 
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where #gc is the observed scattering diameter of 
Sgr A*. If the scattering screen is, in fact, associated 
with the GC (e.g., Dqq < 100 pc), then extragalac- 
tic sources viewed through the screen necessarily have 
angular diameters that exceed 1 arc min at 1 GHz. 

The screen geometry we have utilized is highly ide- 
alized. Other, more physical, geometries are also pre- 
sented in Fig. |l|. The scattering screen may surround 
the GC with spherical or cylindrical geometry. Al- 
ternately, the scattering region may be distributed 
throughout the GC region rather than being confined 
to a thin screen. Yusef-Zadeh et al. (1994), for ex- 
ample, interpret anisotropic scattering from Sgr A* 
in terms of H n regions distributed uniformly within 
100 pc of Sgr A*. The scattering from Sgr A* may 
then be effectively screen-like because it will be dom- 
inated by H ii regions farthest from Sgr A* (through 
geometric effects). However, background sources will 
be scattered by all H n regions. In these cases, the 
angular diameters of background sources are at least 
double what we have derived if the scattering is from 
a thin screen that surrounds the GC, but it will be 
even greater for scattering regions that fill the GC. 

Van Langevelde et al. (1992) derived a lower bound 
on the GC-screen distance based on the fact that the 
free-free optical depth must be ^ 1 at 1.6 GHz: 
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In this equation, £q, t\ are the outer and inner 
scales, respectively, of the electron-density fluctua- 
tions; g(y,'l') is the Gaunt factor, and T4 is the tem- 
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Fig. 1. — Geometries for the Galactic Center scat- 
tering region. (Top) A single, thin screen located a 
distance Agc from the Galactic center. This geom- 
etry is used in most of the dicusssion in this paper. 
(Middle) A spherical (or cylindrical) shell surround- 
ing the GC. (Bottom) A scattering region composed 
of "screenlets" distributed around the GC with total 
size ~ A GC . 

can be adjusted to reduce the overall coefficient from 
0.1 to 0.004 (e.g., through an inner scale of 0.01 pc and 
a temperature of 10 5 K). However, even 10 5 K may 
be too hot for electron de nsity variations to be sus- 
tained (e.g., Spanglcr 1991 ). A scattering screen thin- 
ner than 0.01 pc is also improbable because s tellar- 
wind shocks tend to be at least this thickness (Kahn 

The outer scale 



Breitschwcrdt 199C) 



which is 

probably comparable to the screen thickness, can- 
not be reduced much below 0.01 pc. Other evidence 
also indicates that the inner scale is not less than 
100 km. This includes the observed scaling of scat- 
tering diameters as v~ 2 for Sgr A* and other heav- 



ily scattered sources such as Cyg X-3 (Molnar et 



perature in units ot ID K. The various parameters 



al. 1995; Wilkinson, Narayan, & Spencer 1994) and 
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NGC6334B flMoran et al. 199C| ). We therefore take 
as a hard lower bound, A GG > 33 pc. This suggests 
that a strong upper bound on the sizes of extragalac- 
tic sources is 9 xga \ < 11' for a screen- like geometry 
that wraps around Sgr A*. 

It is possible to determine the screen distance from 
the Galactic center by either measuring the diameters 
of extragalactic sources or, in the case where severe 
scattering diminishes our ability to detect sources, 
use the deficiency of sources to constrain the screen 
location. Attempts to observe extragalactic sources 
near Sgr A* in fact show a deficiency of sources 
within 0.5° of Sgr A* flLazio ct al. 1996|). A likeli- 



ho od analysis shows that this deficiency is consistent 



in the free electron density that are physically present 
(cf. |Cordes fe Lazio 1991 ). This is consistent with 
the observed v~ 2 scaling of the angular diameter of 
Sgr A* . The geometric factor is / — > 1 if the screen is 
midway along the line of sight. But for screens very 
near the GC, / — ► x^ 1 3> 1. Therefore, pulsars at 
the same location as the GC will show at least 6.3 s 
of pulse broadening at 1 GHz^|. The pulse broaden- 
ing may be significantly larger, perhaps as much as 
200 times larger, because the scattering screen may 
be only 33-100 pc from the GC. The minimal scat- 
tering time of 6.3 s may be compared, at 1 GHz, to 
the pulse broadening o f the most heavily scattered 
pulsar, PSR B1849-00 ([Frail fc Clifton 1989j |Clifton 



wit h extragalactic source sizes larger than 0.5' and et al. 1987[ ), which is about 0.3 s 



A GC <; 200 pc (|Lazio ct al. 1996[) 



3. PULSE BROADENING 

Pulse broadening is a diffraction phenomenon but 
may be treated with stochastic ray tracing through 
extended media (Williamson 1972, 1975). It has been 
observed for many pulsars and used to study the dis- 
tribution of ionized microturbulence in the Galaxy 
( jCordes et al. 1991| ). 

The pulse broadening due to the screen responsible 
for the scattering of Sgr A* has an e _1 time scale 



(3) 



where the screen's location along the line of sight is 
represented by the geometric factor^] 



f{x) = x-\l-x). 



(4) 



The pulse broadening time for fiducial values of the 
distance (8.5 kpc) and scattering diameter (1.3") at 
a frequency of 1 GHz is 
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In Eq. U we have adopted a frequency scaling cx v~ A 
rather than the often encountered v scaling be- 
cause, in the extremely strong scattering limit, the 
scattering is dominated by the smallest irregularities 



^"Pulse broadening is often expressed in terms of the screen scattering 
angle 6 S rather than the observed angle <9 Q . Using 9 S , the equivalent 
geometric factor is x(l — x) , which maximizes at x — 1/2. 



Pulsars beyond the GC (but still behind the scat- 
tering screen) will show even larger scattering. For a 
pulsar distance D > -Dec — A G c, the pulsar-screen 
distance is A = D — £> GG + A GG and the pulse broad- 
ening from the screen is 



tgc(D) ~ r GC (D GC ) (^p) 



'GC , 



(6) 



As a function of distance from the Sun, pulse broad- 
ening increases slowly and according to the TC model, 
which possesses components that grow stronger in the 
inner Galaxy. Then, just beyond the location of the 
GC scattering screen, pulse broadening increases dra- 
matically and continues to increase. To combine the 
TC model and the GC screen component, we write 
the net pulse broadening as 



ttc, 



D < Dqc - A GC ; 



(ttc+4c) 1/2 , D>D GC -A 



(7) 
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Combining the TC and GC-screen scattering times 
is ad hoc in form but is sufficiently accurate for our 
purposes here because the GC component is much 
larger than the TC contribution. 

Figure || shows the pulse broadening at two fre- 
quencies (1.4 and 10 GHz) for a range of GC-screen 
distances, A GC = 0.05, 0.1, 0.2, 1.0, and 4.25 kpc. For 
pulsars beyond the GC, the pulse broadening asymp- 
totes to r GC (L» GC )(L» GC /A GC ) ~ 10 5 s at 1.4 GHz 
and 18 s at 10 GHz. 



An early analysis (Davics, Walsh & Booth 1976) estimated 10 s of pulse 
broadening at 1 GHz while implicitly assuming the scattering region to 
be midway between us and the GC. 
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Fig. 2. — Pulse broadening is plotted against dis- 
tance for 5 separate values of the GC-screen dis- 
tance, Agc: 0.05 kpc (thickest line), 0.1, 0.2, 1.0, and 
4.25 kpc (thinnest line). The left-hand scale applies 
to 1.4 GHz, the right-hand scale to 10 GHz. Broaden- 
ing at other frequencies may be estimated using the 
assumed z^~ 4 scaling. 

4. DETECTION OF SCATTERED PULSARS 
IN PERIODICITY SEARCHES 

Pulse broadening decreases the number of har- 
monics that exceed a predetermined threshold in the 
power spectrum of the intensity, thus reducing the 
sensitivity of a pulsar search. Consider a train of 
pulses with period P, average pulse area A , duty 
cycle e, and pulse width (FWHM) W = eP. The dis- 
crete Fourier transform of the pulse train is a series 
of spikes at frequencies t/P, I = 0, 1, . . . each having 
an amplitude, 



A DFT (e) = A ~g(ee), 



(8) 



where pulses have a generic shape g(4>) in pulse phase 
</>, whose continuous Fourier transform is g. 



We define the intrinsic pulsed fraction of the pulsar 
flux as the ratio of the fundamental frequency and 
zero frequency ("DC") amplitudes: 
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where the third equality is for gaussian-shaped pulses 
[i.e., g{4>) — exp(— 41n20 2 )]. For most pulsars, e < 
0.1 implying rjp ~ 1. 

Broadening increases the pulse width to W c g ~ 
(W 2 + t 2 ) 1 / 2 and, hence, the duty cycle to e off = 
W e g/P. The pulsed fraction becomes 
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Figure || shows the pulsed fraction plotted against 
frequency for five different pulse periods. We have 
assumed that the GC screen is near Sgr A* (Agc = 
50 pc). 

Most pulsars have radio spectra that are power- 
laws in form, some objects displaying one or more 
break points separating power-laws of different slope 
( Lorimcr ct al. 1995| ). Letting the spectrum be 

— fs) 

S v cx v we maximize the pulsed flux density rfp S v 
against frequency to solve for the corresponding opti- 
mal pulse broadening, r max = (aln2/2) 1 / 2 P/7r. The 
frequency at which the pulsed flux density is maxi- 
mized can be found from Eq. |^ as 
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= 2.41GHz 
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where Ao.i = Acc/O.lkpc. Note that the optimal 
frequency is a function of pulse period. At the optimal 
frequency, the pulsed fraction is r]p S ^ — e~ a ^ 8 and the 
pulse broadening yields an effective duty cycle 



e c ff = (e 2 +aln2/27r 2 



,1/2 



0.19Va. 



(12) 



The approximate equality is for the case where e <C 
0.2 and suggests that, for such pulsars, the search 
algorithm should search for no more than e 



harmonics in the Fourier transform. Figure W shows 
Vp Sv plotted against frequency for several pulse pe- 
riods, P, and spectral indices, a. 

There is a barrier of indetectability at low fre- 
quencies. It is obvious that only very luminous, 
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Fig. 3. — Plot of the pulsed fraction rjp against fre- 
quency v for different pulse periods, P = 5 s (thickest 
line), 1 s, 0.3 s, 0.1 s, 30 ms, & 5 ms (thinnest line). 
The scattering screen is assumed to be close to the 
GC, Aqc = 50 pc. If the screen is more distant from 
the GC than 50 pc, the plotted curves move toward 
the left. 

long period pulsars with shallow spectra are candi- 
dates for detection. The most luminous pulsars are 
those that show "core" emission ( Rankin 199C ). A 
good example is B1933+16 with P = 0.36 s, dis- 
tance D — 7.9 kpc, and period-averaged flux den- 
sities S400 = 242 mJy and 6*1400 = 42 mjy at 400 and 



1400 MHz, respectively (Taylor et al. 1993). If placed 



at the GC, the pulsed fraction of this pulsar > e 
for v > 2.4/ 1 / 4 GHz. With a = 1.4 and / = 170 (i.e., 
Agc = 50 pc), we find that ^ max « 10.8 GHz. At this 
frequency, the pulsed flux of B1933+16 is 1.8 mJy if 
placed at the GC. Pulsar surveys can reach this flux 
density through use of large bandwidths (e.g., 1 GHz). 
For pulsars in general, we find that the optimal range 
of search frequencies is 5 GHz <; v max <; 18 GHz for 
0.5 £ a & 2 and 0.1 £ P £ 5 s. 



Fig. 4. — The pulsed flux, taking scattering into 
account, as a function of frequency. Pulsar spectra 
of the form v~ a have been normalized so that the 
flux at 1 GHz is the same for all a. The intrinsic 
pulse duty cycle is e = 0.05. Solid Lines: P = 1 s for 
a = 0.5 (heaviest line), 1, 1.5, 2, 2.5 (thinnest line). 
Dotted Line: P = 5 s and a — 2.5. Dashed Line: 
P = 30 ms and a = 0.5. 

The frequency of maximum, pulsed flux, as we have 
defined it, does not correspond to the frequency where 
S/N is maximized in a Fourier transform search. Max- 
imizing S/N depends on the number of harmonics 
detectable in the DFT, and on the variation of sys- 
tem temperature with frequency. At the optimal fre- 
quency, where approximately ~ 5 DFT spikes 
should be detectable, the S/N can be increased by a 
factor ~ \fh. 

We conclude this section with an example search 
program for finding GC pulsars in periodicity searches 
(see Table [l]). We consider searches at 5, 8, 10.8 
and 15 GHz on a telescope with sensitivity equal to 
that for the National Radio Astronomy Observatory's 
Green Bank Telescope, now under construction. The 
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minimum flux density for a single harmonic in the 
DFT is 



130 M Jy, 



(13) 



where T sys is the system temperature, assumed to be 
50 K, G ~ 2 K/Jy is the telescope gain, the receiver 
bandwidth is Az/ = 1 GHz, and the total observing 
time is T =1 hr. The signal-to-noise threshold for de- 
tection of a signal is assumed to be t\t — 10. The 1- 
GHz bandwidth is substantially larger than that used 
in most surveys but is quite feasible at high frequen- 
cies. Moreover, only coarse channelization is required 
to combat dispersion smearing. For dispersion mea- 
sures as large as 2000 pc cm -3 , only about 16 channels 
are needed. The results in Table |l| illustrate an analy- 
sis of the fundamental only; as discussed above, anal- 
ysis of higher harmonics could improve the quoted 
sensitivities by a factor of ~ 2. 

Two different measures of a search's ability to de- 
tect pulsars are illustrated in Table [j]. The first 
is its ability to detect B1933+16, if it were placed 
in the GC. Entries in the Table suggest that large 
S/N can be achieved at frequencies greater than 8 
GHz. Detection at 5 GHz is hopeless. We also 
assess the fraction of the total GC pulsar popula- 
tion that a search could detect, which we determine 
from the luminosity function for radio pulsars. In 
the solar neighborhood, for observation frequencies 
near 400 MHz, the pulsar luminosity function has the 
form dN/d\ogL oc with /? « 1 (Lyne, Manch- 
ester, & Taylor 1985; [Manchester fc Taylor 1977[). 



The fraction of sources with luminosity greater than 
L is /l(> L) w (1 — L/Li)L /L, where Lq,L\ are 
the minimum and maximum intrinsic luminosities of 



a p ulsar. In th e most recent catalog of pulsars (Tay- 



lor ct al. 1993 ), only one pulsar has a luminosity less 
than 1 mJy kpc 2 , so we take Lq = 1 mJy kpc 2 . We 
use a spectral index of 1.4 to scale in frequency. For 
S m in given by Eq. the implied minimum luminos- 
ity at the particular observation frequency is ~ 11 
mJy kpc 2 . Scaling from the frequencies in the Table 
then yields the equivalent L at 400 MHz and the frac- 
tion of objects that exceed this minimum. We also 
tabulate the period, P, at which the tabulated fre- 
quency is, in fact, the optimal frequency, as given by 
Eq.[ll]. As can be seen, this period exceeds that of any 
known radio pulsar at 5 GHz, but decreases progres- 
sively to ~ 0.1 s at 15 GHz. The fraction of objects 
that exceed the minimum luminosity also decreases in 
going from 5 to 15 GHz. A key disadvantage to us- 



ing higher frequencies is that, all other things being 
equal, the telescope time needed to search a region of 
fixed angular size will scale as the solid angle of the 
telescope beam, oc v~ 2 . This amounts to nearly an 
order of magnitude increase in time needed in going 
from 5 to 15 GHz. 

The results in Table [l] are only illustrative. Pulsars 
show a wide range of spectral indices and there are 
possible correlations between period and luminosity 
that we have ignored. Also, we have been assuming 
a constant pulse-broadening value for pulsars in the 
GC. In actuality, the spatial distribution of pulsars in 
the GC will create a range of pulse broadening. The 
extreme values of pulse broadening we have estimated 
above will be obtained if the angular distribution of 
pulsars is concentrated toward Sgr A*, yet are still 
quite far from the scattering screen. However, pul- 
sars spread throughout the region will show a wide 
range of broadening. Objects close to the screen will 
show small scattering, asymptotic to the TC value 
as A psr — > 0. Periodicity searches at different (high) 
frequencies will therefore sample the population at 
different A psr . 

5. SEARCHES FOR PULSARS IN IMAG- 
ING SURVEYS FOR POINT SOURCES 

Scattering preserves the total flux density while 
attenuating the pulsed flux. Pulsars may therefore 
be searched for as objects with pulsar-like spectra 
and polarization in aperture synthesis surveys of the 
GC. Angular broadening dilutes the surface bright- 
ness, but we show here that this dilution denigrates a 
synthesis survey far less than pulse broadening does 
a periodicity search. Indeed Sgr A*, the OH/ IR stars 
(sQ), and the GC transient ( Zhao et al. 1992 ) demon- 



strate that heavily scattered, yet sufficiently lumi- 
nous, sources in the GC can be detected with aperture 
synthesis instruments. 

Table ^ illustrates detection of pulsars in an imag- 
ing survey with the VLA. We assume observations are 
conducted with a total integration time of 1 hr and 
total bandwidth At/ = 50 MHz at 1.4 and 5 GHz and 
Av = 3 MHz at 0.33 GHz. As before, we consider 
the detection of both B1933+16 and a fraction of the 
total pulsar po pulation Re cent observations at 

0.33 GHz (e.g., Frail et al. 199?;) have demonstrated 
the capability of making low-frequency observations 
with the VLA which approach the thermal noise limit 
in the A-configuration. Within Sgr A West, i.e., 
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Table 1 

Pulsar Detection in Periodicity Searches 
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Note. — (1) Observation frequency; (2) Pulse-broadening time for Aqc = 50 pc; (3) S sys — T sys /G, where G is 
the telescope sensitivity in K Jy -1 ; (4) Minimum detectable flux of fundamental in FFT; (5) Pulsed flux density 
of PSR B1933+16; (6) Signal-to-noise ratio for a search of PSR B1933+16, if it were at the Galactic center; (7) 
Period for which the frequency is the optimum search frequency. (8) Fraction of total pulsar population likely to 
be detectable, assuming nominal spectral index of —1.4 



within 5 arc min of Sgr A*, strong thermal absorp- 
tion is seen at 0.33 GHz. This absorping gas strongly 
attenuates Sgr A* below 1 GHz ( Anantharamaiah ct 
al " 



1991; Pedlar et al. 



and may contribute to 
the extreme scattering of Sgr A* at higher frequen- 
cies. However, heavily scattered OH masers are seen 
up to 25 arc min from Sgr A* (Frail et al. 1994). 
Over this larger region, the frcc-frcc optical depth is 
only ~ 1. Therefore, absorption may attenuate source 
fluxes but should not render the entire scattering re- 
gion opaque. 

For the Table, we have assumed that pulsar scat- 
tering diameters are the same as that of Sgr A* . From 
considerations in §|[ the scattering diameter of a pul- 
sar at distance D > Dgc ~ Aqc is 
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Pulsars closer to the screen than Sgr A* will show less 
angular broadening than Sgr A* while a pulsar only 
Aqc further than Sgr A* will be twice the diameter. 
An object 0.5 kpc beyond Sgr A* will be 10 times 
the diameter (for Aqc = 50 pc), or 13"^^. The 
angular diameters of pulsars in the vicinity of Sgr A* 
are therefore well matched to the synthesized beams 
of the VLA in its four configurations, though only the 
larger configurations filter out the intense extended 



emission sufficiently. Note that the pulsar B1933+16 
is easily detectable if placed in the GC at 1.4 as well 
as 5 GHz. Moreover, about 1.2% of all pulsars will be 
detectable at 1.4 GHz (for an assumed spectral index 
of 1.4), about 10 times larger than the fraction found 
in a periodicity search. 

6. GALACTIC CENTER PULSAR POPU- 
LATIONS 

The preceding sections have assumed the existence 
of pulsars in the GC. In this section we assess how 
many pulsars are likely to exist in the GC and use the 
results of §§|| and || to indicate how many of these are 
in fact detectable. The possibility of a large number 
of stellar remnants at the GC has been discussed by a 
number of other authors, though their focus has been 
on the central few parsecs (e.g., Saha, Bicknell, & Mc- 



Gregor 1996) and on either black holes (Morris 1993) 
or white dwarfs ( Haller et al. 199"6| ). Our focus will 
be on neutron stars exclusively and, based on the size 
of the scattering region, will be on a larger size scale 
than considered previously. We estimate the number 
of GC pulsars by considering both the integrated star 
formation history of the GC and the effects of a recent 
star burst. Finally, we also discuss the implications of 
the detection of GC pulsars on our understanding of 
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Table 2 

Pulsar Detection in an Imaging Survey 



B1933+16 



V 0GG /min I (S/N) f L 

(GHz) (") (mJy/beam) (mJy/beam) 

(1) (2) (3) (4) (5) (6) 



0.3 11 2.1 a 40 19 0.007 

1.4 0.62 0.17 33 194 0.012 

5 0.05 0.13 6.7 52 0.002 



Note. — (1) Observation frequency; (2) Scattering angle for Galactic center source (i.e., Sgr A*); (3) Minimum 
detectable brightness; (4) Brightness of PSR B1933+16, if it were at the Galactic center; (6) Signal-to-noise ratio 
for a search of PSR B1933+16, if it were at the Galactic center; (7) Fraction of total pulsar population likely to be 
detectable, assuming nominal spectral index of —1.4 

a Assumes bandwidth at 0.33 GHz of Av = 3 MHz. 



the dynamics and star formation history of the GC. 

We first consider the number of pulsars likely to 
have been formed if the star formation in the GC has 
been occurring at a constant rate over the Galaxy's 
history. In actuality, factors such as episodic star 
bursts, an initial mass function (IMF) different from 
that in the disk, and the shape of the central po- 
tential and its effects on gas dynamics may all have 
contributed to a variable star formation rate. The as- 
sumption of a constant star formation rate is nonethe- 
less useful for illustration and comparison with the 
starburst model. 

6.1. Steady Star Formation 

We extrapolate from the current mass of the GC, 
for an assumed IMF, in order to estimate the num- 
ber of pulsars likely to have formed within the GC. 
Various methods, summarized by Genzel, Hollenbach, 

6 Townes (1994), constrain the mass enclosed within 
the central 100 pc to be in the range ~ 4 x 10 8 - 
10 9 Mq. An estimate of the enclosed mass derived 
from the stellar luminosity distribution, and assum- 
ing a constant mass-to-light ratio, is approximately 

7 x 10 s Mq, which we adopt as a nominal value for 
the mass interior to 100 pc. Most o f this mass is in 
the form of stars ( Genzel et al. 1994 ). Because we do 
not know the form of the IMF in the GC, we consid- 



ered different IMFs, one a ppropriate for the Galactic 
disk ( |Miller fc Scalo 1979]) and one that accounts for 
a stronger weighting of the IMF toward high-mass 
stars during starbursts ( Rieke et al. 198C ). For these 
IMFs, and different assumed ranges of masses for 
main-sequence stars that result in neutron stars, we 
find that 4-12 x 10 7 neutron stars have formed in the 
GC. 

Implicit in this extrapolation is the assumption 
that star formation and evolution in the GC is sim- 
ilar to that in the disk. Morris (1993) has discussed 
the conditions in the central 1-10 pc and how they 
could both inhibit massive star formation and alter 
the evolution of those stars which do form. Meynet et 
al. (1994) have shown that massive, high-metallicity 
stars may lose a substantial fraction of their mass dur- 
ing their post-main-sequence evolution and they sug- 
gest the remnants of such stars would be white dwarfs 
rather than neutron stars (and black holes). Even if 
such processes are operative, however, the presence 
of some massive stars is clearly indicated. There are 
a number of H n regions within 15' of Sgr A* which 
appear to be powered by OB a nd/or WR stars (e.g ., 
Figer, McLean, & Morris 1996; [Poglitsch et al. 1996| ); 
observations of H2 (v = 1 — * 0) S(l) emission in 
the inner 2° imply excitation by a UV radiation field 
(Pak, Jaffe, & Keller 1996); and star counts within 
the inner 2° reveal a number of stars with dered- 
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dened absolute magnitudes of Mk ^ —9, suggestive 
of stars with M ^ 5 Mq (Catch pole, Whitelock, 
& Glass 1990; |Bertelli et al. 1994| ). Bearing these 
caveats in mind, we shall use our estimate of 10 7 -10 8 
neutron stars in the GC. 

The actual number of active and detectable pulsars 
in the GC will be substantially smaller. Firstly, the 
mean space velocity of pulsars is 500 km s _1 (Lyne & 



Loiimer 1994). A considerable fraction, ;> 0.25, is not 
bound to the Galaxy. Only those pulsars forming the 
lowest-velocity tail, < 100 km s _1 or about 16% of the 
objects, are bound to the central 100 pc. Secondly, 
the lifetime for pulsar radio emission is ~ 10 7 yr for 
strong field objects with surface fields ~ 10 12 Gauss. 
If the Galactic bulge is ~ 10 10 yr old, only a fraction ~ 
10~ 3 of the neutron stars formed in the central 100 pc 
will still be active pulsars. Hence, of the ~ 10 7 -10 8 
neutron stars formed in the GC over the lifetime of the 
Galaxy, ~ 10 3 -10 4 are likely to be active pulsars still 
within the central 100 pc. Of thes e, only a fr a ction 
/b ^0.2 will be beamed toward us (Biggs 1990|; Lync 
& Manchester 1988| ; [Marayan fe Vivekanand 1983|) . If 
about 1% of these can be detected, as in an imaging 
survey using the VLA, then only a few to ~ 20 objects 
are expected to be found. A periodicity search that 
finds only the 0.1% most luminous pulsars is unlikely 
to find any pulsars from this population. 

One effect which would increase this number is the 
production of recycled pulsars in binaries resulting 
from tidal captures and stellar collisions. The cen- 
tral stellar density may exceed 10 pc~ 3 , sufficiently 
high that stellar encounters should be frequent and 
that processes for recycling pulsars should be oper- 
ative. Fabian, Pringle, & Rees (1975) estimate the 
tidal capture rate (see also Hut ct al. 1992). If the 
number of neutron stars formed in the GC has been 
~ 10 7 -10 8 , then ~ 10 3 tidal captures or collisions will 
have occurred over the Galaxy's history. Not all of 
the binaries so formed will produce a recycled pulsar 
and primordial binary systems could also contribute 
to the population of recycled pulsars. Regardless of 
formation mechanism, the detection of a large num- 
ber o f soft, compact X-ray sou rces toward the GC 
(e.g., Predehl fc Trumpcr 1994 ) supports the possi- 



bility of an enhanced number, relative to the disk, 
of neutron-star binary systems in the GC. Although 
some of these X-ray sources are probably black-hole 
systems, the fraction of millisecond pulsars in the to- 
tal pulsar population in the GC could nevertheless be 
enhanced substantially relative to a value ~ 0.5 in the 



disk ( Bailes fc Lorimcr 1995| ; Lorimcr ct al. 1993| ). 

Unfortunately, the fastest spun-up pulsars will not 
be detectable in a periodicity search, though they may 
appear in an imaging search. Modestly spun-up pul- 
sars with lifetimes 10-100 times longer than those of 
most pulsars may be the best targets. In this re- 
spect it is interesting that 15% of the spun-up pulsars 
in a recent listing had periods greater than 100 ms 
(Kulkarni 1995), though none of these objects would 
themselves be detectable if placed at the GC. 



6.2. Star Bursts 

The discussion above assumes a constant star for- 
mation rate within the GC. A recent starburst has 



been invoked (e.g., Hartmann 1995; Gcnzcl et al. 1994; 
Sofue 1994; but see also Morris 1993) to explain one 
or more features of the GC. Estimates for the age 
of the burst are 10 6 -10 7 Myr with the total num- 
ber of neutron stars formed being 10 2 -10 6 . Even if 
a burst per se has not occurred, several lines of evi- 
dence point to an enhanced star formation rate, rel- 
ative to the disk rate, during the recent past in the 
GC. These include the density of supernova remnants 
within the central 10° which is double that found in 



the disk ( Gray 1994 ) , and the presence of a large num- 
ber of shell and arc structures in molec ular gas, also 
suggestive of recent supernova activity ( Hasegawa et 
al. 1996|) . 

The estimated range of ages for the most recent 
starburst is sufficiently small that all of the pulsars 
formed during the burst should still be active. High 
velocity pulsars will escape the GC in a time short 
compared to the burst. Indeed, Hartmann (1995) ad- 
vocates searches for these escaped pulsars as a means 
for diagnosing the recent star formation history in the 
GC. Depending upon the number of neutron stars cre- 
ated in the burst, the total number of pulsars in the 
GC from the low-velocity tail of the pulsar velocity 
distribution could number 10-10 5 . That is, the num- 
ber of active pulsars from the burst population could 
exceed that from a constant star formation rate over 
the history of the GC. 

The number of detectable pulsars will be smaller 
than the total number active and within the GC 
due to selection effects. In §§|] and || we estimated 
the fraction of pulsars sufficiently luminous to be de- 
tected, /l, based on the local pulsar luminosity func- 
tion. This fraction is 10~ 2 -10 -5 , depending upon the 
search method and frequency. Taking into account 
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the non-isotropic nature of pulsar emission, the re- 
sulting fraction of detectable pulsars is ~ 10~ 3 -10~ 6 . 
Thus, the number of active, detectable pulsars from 
the central 100 pc of the Galaxy could be <^ 1 to as 
many as 100. A constant star formation rate over the 
GC's history will likely contribute <; 10 pulsars to the 
total number detectable. 

6.3. Pulsars as Probes of the Galactic Center 

We conclude this section with a brief discussion of 
the utility of pulsars for constraining the dynamics 
and star formation history of the GC. The number 
of pulsars detected, together with a careful account- 
ing for selection effects, can be used to constrain the 
past star formation rate of the GC. If any pulsars 
can be detected as sources of pulsed emission, the 
distribution of spin periods; their period derivatives, 
P; and dispersion measures, DM, and pulse broad- 
ening should also be measurable. The distribution 
of periods can be used in estimating the past star 
formation rate (e.g., Lorimer et al. 1993 ). Measure- 
ment of P has the potential to constrain the GC 
gravitational potential, much like the determination 
of a negative P for B2127+11 has set limits on its 
acceleration from the potential of the globular clus- 
ter Ml 5 flWolszczan et al. 1989] ) . As noted in §0, 
the most recent model for the Galactic distribution 
of ionized gas has few constraints toward the inner 
Galaxy. Even a few DM and pulse broadening mea- 
surements would prove useful for improving substan- 
tially our knowledge of conditions there. Moreover, 
gas velocities, particularly near Sgr A*, are large, in 
some cases exceeding 100 km s _1 (Schwarz, Bregman, 
& van Gorkom 1989). DM monitoring programs can 
expect to detect changes in the DM on time scales of 
a month, which in turn could be used to probe the 
AU-scale structure of ionized gas in the GC. Even if 
pulsars can only be detected in an imaging survey, 
their angular diameters would still provide both con- 
straints on the distribution and a probe of the small- 
scale structure of ionized gas in the GC. 

Finally, the clock-like nature of pulsars has prompted 
speculation about their utility for studying the con- 
stituents and mass distribution in the GC. Wex, Gil, 
& Sendyk (1996) and, earlier, Paczynski & Trim- 
ble (1979) have discussed the microlensing and gravi- 
tational time delay of GC pulsars. The time delay, in 
particular, can be used to constrain the mass of any 
central object. Such determinations, however, will be 
possible at only the highest frequencies where pulsed 



emission is manifested. 

7. CONCLUSIONS AND RECOMMENDA- 
TIONS 

The extreme level of scattering seen toward the 
Galactic center (1.3" of broadening for Sgr A* at 
1 GHz) produces considerable pulse broadening. The 
amount of pulse broadening depends upon the Galac- 
tic center-scattering screen distance, Aqcj but at 
minimum, the pulse broadening is 6.3 seconds. If the 
screen is closer to the Galactic center, the pulse broad- 
ening is increased. Radio pulsars within about 15'-1° 
of Sgr A* and behind the GC scattering screen will 
only be detected by periodicity searches conducted at 
frequencies much higher than used heretofore or by 
imaging surveys. 

Periodicity searches at 5 GHz will be capable of 
detecting pulsars with canonical spin periods ~ 1 s, 
though with much reduced sensitivity if Aqc "C Dqq ■ 
At 10 GHz, pulse broadening of millisecond pulsars is 
a significant fraction of a period even if Aqc ~ Dgg . 
If Agc ~ 100 pc, pulsars with periods shorter than 
0.5 s at 5 GHz and 34 ms at 10 GHz will be inacces- 
sible. To reach pulsars with the Crab pulsar's period 
(33 ms) such that a 10% duty cycle pulse may be seen 
without broadening, the observing frequency must be 
at least 18 GHz. While the Crab pulsar itself is unde- 
tectable at this frequen cy, other young pul sars with 
shallower spectra (e.g., Lorimer et al. 1995 ) may be 
strong enough. Such objects are rare, so we conclude 
that imaging surveys at fairly low frequencies (e.g., 
1.4 GHz) are much more likely to be successful in 
finding Galactic center pulsars. 

The optimal frequency for conducting any period- 
icity searches is that which maximizes the pulsed flux 
density and is given by v max = 7.3 GHz (y^aAo.iP) 1 ^ 4 
(Eq. [ll]). The effective duty cycle of a scattered pul- 
sar at the optimal frequency is ~ and suggests 
that search algorithms need not consider more than 



c cir 



5 harmonics in the DFT. 



The angular broadening for a pulsar at the GC 
will be similar to that of Sgr A* ; the broadening for a 
pulsar more distant than Sgr A* but seen through the 
same scattering screen is given by Eq. |l4|. These an- 
gular sizes are well matched to the synthesized beams 
of the VLA in its four configurations. 

Finally, the number of pulsars detectable in the 
Galactic center depends upon the star formation his- 
tory of the Galactic center. The largest number of 
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detectable pulsars will occur if the GC has undergone 
a recent episode (<; 10 7 yr) of massive star forma- 
tion. The number of active, detectable pulsars in the 
GC may be as many as 100. Smaller starbursts and 
the star formation over the Galaxy's history can con- 
tribute another 1-10 pulsars. 

We thank Z. Arzoumanian, D. Chernoff, D. Hart- 
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